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Using multiple stochastic integrals and the Malliavin calculus, 
we analyze the asymptotic behavior of quadratic variations for a spe- 
cific non-Gaussian self-similar process, the Rosenblatt process. We 
apply our results to the design of strongly consistent statistical esti- 
mators for the self-similarity parameter H . Although, in the case of 
the Rosenblatt process, our estimator has non-Gaussian asymptotics 
for all H > 1/2, we show the remarkable fact that the process's data 
at time 1 can be used to construct a distinct, compensated estimator 
with Gaussian asymptotics for H £ (1/2,2/3). 

1. Introduction. 

1.1. Context and motivation. A self-similar process is a stochastic pro- 
cess such that any part of its trajectory is invariant under time scaling. 
Self-similar processes are of considerable interest in practice in modeling 
various phenomena, including internet traffic (see, e.g., [32]), hydrology (see, 
e.g., [13]) or economics (see, e.g., [12, 31]). In various applications, empirical 
data also shows strong correlation of observations, indicating the presence, in 
addition to self-similarity, of long-range dependence. We refer to the mono- 
graphs [7] and [25] for various properties of, and fields of application for, 
such processes. 

The motivation for this work is to examine non-Gaussian self-similar pro- 
cesses using tools from stochastic analysis. We will focus our attention on a 
special process of this type, the so-called Rosenblatt process. This belongs to 
a class of self-similar processes which also exhibit long-range dependence and 
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which appear as hmits in the so-cahed noncentral Umit theorem: the class 
of Hermite processes. We study the behavior of the quadratic variations for 
the Rosenblatt process Z, which is related to recent results by [14, 16, 17], 
and we apply the results to the study of estimators for the self-similarity 
parameter of Z. Recently, results on variations or weighted quadratic vari- 
ations of fractional Brownian motion were obtained in [14, 16, 17], among 
others. The Hermite processes were introduced by Taqqu (see [27] and [28]) 
and by Dobrushin and Major (see [5]). The Hermite process of order q>l 
can be written, for every t>0, as 



Zj,{t)=c{H,q) 

(1.1) 



[Jo Vti / . 



dW{yi)- --dW ivg), 



where c{H, q) is an explicit positive constant depending on q and H , and such 
that E(Z|f (1)^) = 1, = max(j;, 0), the self-similarity (Hurst) parameter H 
belongs to the interval (i, 1) and the above integral is a multiple Wiener-Ito 
stochastic integral with respect to a two-sided Brownian motion iW {y))y^^ 
(see [21]). We note that the Hermite processes of order g > 1, which are non- 
Gaussian, have only been defined for H > ^] how to define these processes 
for i/ < i is still an open problem. 

The case g = 1 is the well-known fractional Brownian motion (fBm): this 
is Gaussian. One recognizes that when g = l, (1-1) is the moving average 
representation of fractional Brownian motion. The Rosenblatt process is the 
case q = 2. All Hermite processes share the following basic properties: 

• they exhibit long-range dependence (the long-range covariance decays at 
the rate of the nonsummable power function it?^~'^); 

• they are //-self-similar, in the sense that for any c > 0, {Z'^{ct))t>o and 
{c^ Zjj{t))t>o are equal in distribution; 

• they have stationary increments, that is, the distribution of {Zjj{t + h) — 
Zjj{h))t>o does not depend on /i > 0; 

• they share the same covariance function, 

B[Zjj{t)Zjj{s)] =: R^{t, s) = i(i2^ + s^^ - \t - s\^^), s,t>0, 

so, for every s,t>0, the expected squared increment of the Hermite pro- 
cess is 

(1.2) B[{ZUt)-Zj,{s)f] = \t-sr 

from which it follows by Kolmogorov's continuity criterion, and the fact 
that each L''(r2)-norm of the increment of Zj^ over [s,t] is commensurate 
with its L2(0) -norm, that this process is almost surely Holder continuous 
of any order S < H; 
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• the qth Hermite process lives in the so-cahed g'th Wiener chaos of the 
underlying Wiener process W since it is a qth order Wiener integral. 

The stochastic analysis of fBm has been developed intensively in recent 
years and its applications are numerous. Other Hermite processes are less 
well studied, but are still of interest due to their long-range dependence, 
self-similarity and stationarity of increments. The great popularity of fBm 
in modeling is due to these properties and fBm is preferred over higher order 
Hermite processes because it is a Gaussian process and because its calculus 
is much easier. However, in concrete situations, when empirical data attests 
to the presence of self-similarity and long memory without the Gaussian 
property, one can use a Hermite process living in a higher chaos. 

The Hurst parameter H characterizes all of the important properties of 
a Hermite process, as seen above. Therefore, properly estimating H is of 
the utmost importance. Several statistics have been introduced to this end, 
such as wavelets. A:- variations, variograms, maximum likelihood estimators 
and spectral methods. Information on these various approaches can be found 
in the book by Beran [1]. 

In this paper, we will use variation statistics to estimate H. Let us recall 
the context. Suppose that a process {Xt)te[o,i] is observed at discrete times 
{0, . . . , ^^^^,1} and let a be a "filter" of length / > and p > 1 a fixed 
power; that is, a is an / + 1-dimensional vector a = (oq, ai, . . . , a/) such that 
J2q=Q ttqq'' = OforO<r<p— 1 and J2q=Q aqq^ ^ 0. The ^-variation statistic 
associated to the filter a is then defined as 



Vat (A:, a) 

where, for i G {/,... , A^}, 



1 



N -I 



N-l 

E 

i=i 



|K(WI^ 



q=0 



N 



When X is fBm, these statistics are used to derive strongly consistent es- 
timators for the Hurst parameter and their associated normal convergence 
results. A detailed study can be found in [8] and [11] or, more recently, in 
[4]. The behavior of VN{k,a) is used to derive similar behaviors for the cor- 
responding estimators. The basic result for fBm is that, if p > H + j, then 
the renormalized /c- variation Vi\f{k,a) converges to a standard normal dis- 
tribution. The easiest and most natural case is that of the filter a = {1, — 1}, 
in which case p = I; one then has the restriction if < |. The techniques 
used to prove such convergence in the fBm case in the above references are 
strongly related to the Gaussian property of the observations; they appear 
not to extend to non-Gaussian situations. 
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Our purpose here is to develop new techniques that can be appHed to both 
the fBm case and to other non-Gaussian self-similar processes. Since this is 
the first attempt in such a direction, we keep things as simple as possible: we 
treat the case of the filter a = {1, —1} with a fc-variation order = 2 (quadratic 
variation), but the method can be generalized. As announced above, we fur- 
ther specialize to the simplest non-Gaussian Hermite process, that is, the 
one of order 2, the Rosenblatt process. We now give a short overview of our 
results (a more detailed summary of these facts is given in the next subsec- 
tion). We obtain that, after suitable normalization, the quadratic variation 
statistic of the Rosenblatt process converges to a Rosenblatt random vari- 
able with the same self-similarity order; in fact, this random variable is the 
observed value of the original Rosenblatt process at time 1 and the con- 
vergence occurs in the mean square. More precisely, the quadratic variation 
statistic can be decomposed into the sum of two terms: a term in the fourth 
Wiener chaos (i.e., an iterated integral of order 4 with respect to the Wiener 
process) and a term in the second Wiener chaos. The fourth Wiener chaos 
term is well behaved, in the sense that it has a Gaussian limit in distri- 
bution, but the second Wiener chaos term is ill behaved, in the sense that 
its asymptotics are non-Gaussian and are, in fact, Rosenblatt-distributed. 
This term, being of a higher order than the well-behaved one, is responsible 
for the asymptotics of the entire statistic. But, since its convergence occurs 
in the mean-square and the limit is observed, we can construct an adjusted 
variation by subtracting the contribution of the ill-behaved term. We find an 
estimator for the self-similarity parameter of the Rosenblatt process, based 
on observed data, whose asymptotic distribution is normal. 

Our main tools are the Malliavin calculus, the Wiener-Ito chaos expan- 
sions and recent results on the convergence of multiple stochastic integrals 
proved in [10, 22, 23] and [24]. The key point is the following: if the observed 
process X lives in some Wiener chaos of finite order, then the statistic Va^ 
can be decomposed, using product formulas and Wiener chaos calculus, into 
a finite sum of multiple integrals. One can then attempt to apply the crite- 
ria in [22] to study the convergence in law of such sequences and to derive 
asymptotic normality results (or to demonstrate the lack thereof) on the 
estimators for the Hurst parameter of the observed process. The criteria in 
[22] are necessary and sufficient conditions for convergence to the Gaussian 
law; in some instances, these criteria fail (e.g., the fBm case with H > 3/4), 
in which proof of nonnormal convergence "by hand," working directly 

with the chaoses, can be employed. It is the basic Wiener chaos calculus that 
makes this possible. 

1.2. Summary of results. We now summarize the main results of this 
paper in some detail. As stated above, we use quadratic variation with a = 
{1,-1}. We consider the two following processes, observed at the discrete 
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times {i/N}fLQ: the fflm process X = B and the Rosenblatt process X = Z. 
In either case, the standardized quadratic variation and the Hurst parameter 
estimator are given, respectively, by 

,,3, .4 

(1.4) = HM2.i-l. 1}) := \ - ^l„g|(x(±) - 

We choose to use the normalization ^ in the definition of Vat (as, e.g., in [4]), 
although, in the literature, it sometimes does not appear. The -ff-dependent 
constants cj^h (et al.) referred to below are defined explicitly in (3.2), (3.6), 
(3.12), (3.14), (3.21) and (3.33). Here, and throughout, L'^{n) denotes the 
set of square-integrable random variables measurable with respect to the 
sigma field generated by W. This sigma- field is the same as that generated 
by B or by Z. The term "Rosenblatt random variable" denotes a random 
variable whose distribution is the same as that of Z(l). 

We first recall the followings facts, relative to fractional Brownian motion: 

1. if X = B and if G (1/2,3/4), then: 



'A^/ci^z/V/v converges in distribution to the standard normal law; 
(b) \/iVlog(A^) ^^A_ (gjv — H) converges in distribution to the standard 
normal law; 

2. if X = ^ and g g (3/4, 1), then: 

(a) ^Jn^~^^^~/c2jiVn converges in L^(r2) to a standard Rosenblatt ran- 
dom variable with parameter Hq = 2H — 1 ; 

(b) N^^^ \og{N)--^^={Hj^ — H) converges in L^(J7) to the same stan- 
dard Rosenblatt random variable; 

3. if X = S and ii = 3/4, then: 



(a) y^A^/(c']^ ^ log A^) Vat converges in distribution to the standard normal 
law; 

(b) ^/W\ogN —7== (if AT (2, a) — H) converges in distribution to the stan- 
dard normal law. 

The convergences for the standardized Vat's in points 1(a) and 2(a) 
have been known for some time, in works such as [28] and [9]. Lately, 
even stronger results, which also give error bounds, have been proven. 
We refer to [19] for the one-dimensional case and ii G (0, f), [2] for then 
one-dimensional case and ii G [f , 1) and to [20] for the multidimensional 
case and H G (0, |). 

In this paper, we prove the following results for the Rosenblatt process 
X = Z as — > 00: 
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4. if X = Z and if G (1/2, 1), then with c^^h in (3.12), 

(a) A^^~^V7v(2,a)/(c3,_f/) converges in L^(r2) to the Rosenblatt random 
variable Z{1); 

(b) ^^;^log(A^)(ffAr(2,a) — H) converges in L^(J7) to the same Rosen- 
blatt random variable Z(l); 

5. if X = Z and H G (1/2,2/3), then, with ei^n and fi^n in (3.21) and 
(3.33), 

(a) [VAf(2, a) — )(jT-H Z(l)] converges in distribution to the stan- 
dard normal law; 

(b) - [21og(A^)(g - HN{2,a)) - ^Z{1)] converges in distri- 
bution to the standard normal law. 

Note that Z{1) is the actual observed value of the Rosenblatt process at 
time 1, which is why it is legitimate to include it in a formula for an estima- 
tor. Points 4 and 5 are new results. The subject of variations and statistics 
for the Rosenblatt process has thus far received too narrow a treatment in 
the literature, presumably because standard techniques inherited from the 
noncentral limit theorem (and sometimes based on the Fourier transform for- 
mula for the driving Gaussian process) are difficult to apply (see [3, 5, 28]). 
Our Wiener chaos calculus approach allows us to show that the standardized 
quadratic variation and corresponding estimator both converge to a Rosen- 
blatt random variable in L'^{Q). Here, our method has a crucial advantage: 
we are able to determine which Rosenblatt random variable it converges to: 
it is none other than the observed value Z{1). The fact that we are able to 
prove L^(r2) convergence, not just convergence in distribution, is crucial. In- 
deed, when H < 2/3, subtracting an appropriately normalized version of this 
observed value from the quadratic variation and its associated estimator, we 
prove that asymptotic normality does hold in this case. This unexpected re- 
sult has important consequences for the statistics of the Rosenblatt process 
since it permits the use of standard techniques in parameter estimation and 
testing. 

Our asymptotic normality result for the Rosenblatt process was specifi- 
cally made possible by showing that Vn can be decomposed into two terms: 
a term T4 in the fourth Wiener chaos and a term T2 in the second Wiener 
chaos. While the second- Wiener-chaos term T2 always converges to the 
Rosenblatt random variable Z{1), the fourth chaos term T4 converges to 
a Gaussian random variable for H < 3/4. We conjecture that this asymp- 
totic normality should also occur for Hermite processes of higher order q>3 
and that the threshold H = 3/4 is universal. The threshold H <2/3 in the 
results above comes from the discrepancy that exists between a normalized 
T2 and its observed limit Z{1). If we were to rephrase results 4 and 5 above, 
with T2 instead of Z{1) (which is not a legitimate operation when defining 
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an estimator since T2 is not observed), the threshold would be < 3/4 and 
the constant fi^n would vanish. 

Beyond our basic interest concerning parameter estimation problems, let 
us situate our paper in the context of some recent and interesting works on 
the asymptotic behavior of variations (or weighted variations) for Gaussian 
processes, namely the papers [14, 15, 16, 17] and [26]. These recent papers 
study the behavior of sequences of the type 

|:Mx((.-i)/A.))( i^'^/^>-/ir'>/^"^ i), 

where X is a Gaussian process (fractional Brownian motion in [14, 16] and 
[17], and the solution of the heat equation driven by a space-time white 
noise in [26]) or the iterated Brownian motion in [18] and /i is a regular 
deterministic function. In the fractional Brownian motion case, the behav- 
ior of such sums varies according to the values of the Hurst parameter, 
the limit sometimes being a conditionally Gaussian random variable, some- 
times a deterministic Riemann integral and sometimes a pathwise integral 
with respect to a Hermite process. We believe that our work is the first to 
tackle a non-Gaussian case, that is, when the process X above is a Rosen- 
blatt process. Although we restrict ourselves to the case when h = l, we 
still observe the appearance of interesting limits, depending on the Hurst 
parameter: while, in general, the limit of the suitably normalized sequence 
is a Rosenblatt random variable (with the same Hurst parameter H as the 
data, which poses a slight problem for statistical applications), the adjusted 
variations (i.e., the sequences obtained by subtracting precisely the portion 
responsible for the non-Gaussian convergence) do converge to a Gaussian 
limit for He {1/2,2/3). 

This article is structured as follows. Section 2 presents preliminaries on 
fractional stochastic analysis. Section 3 contains proofs of our results for the 
non-Gaussian Rosenblatt process. Some calculations are recorded as lemmas 
that are proven in the Appendix. Section 4 establishes our parameter esti- 
mation results, which follow almost trivially from the theorems in Section 3. 

2. Preliminaries. Here, we describe the elements from stochastic anal- 
ysis that we will need in the paper. Consider Ti., a real, separable Hilbert 
space and {B{ip),ip € W), an isonormal Gaussian process, that is, a centered 
Gaussian family of random variables such that E(B{(p)B('il')) = 7/^)7^. 

Denote by In the multiple stochastic integral with respect to B (see [21] 
and [30]). This /„ is actually an isometry between the Hilbert space Ti®"" 
(symmetric tensor product) equipped with the scaled norm — * tviS)'"' and 
the Wiener chaos of order n which is defined as the closed linear span of the 
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random variables Hn{B{ip)), where 9jG'H,||(/5||-^ = 1 and Hn is the Hermite 
polynomial of degree n. 

We recall that any square-integrable random variable which is measurable 
with respect to the d-algebra generated by B can be expanded into an 
orthogonal sum of multiple stochastic integrals, 

n>0 

where /„ G 7^®" are (uniquely determined) symmetric functions and /o(/o) = 

nn 

In this paper, we actually use only multiple integrals with respect to the 
standard Wiener process with time horizon [0, 1] and, in this case, we will 
always have TL = L'^ {[0,1]). This notation will be used throughout the paper. 

We will need the general formula for calculating products of Wiener chaos 
integrals of any orders, p and q, for any symmetric integrands / S TC-^p and 
g £ n®"; it is 

pAg / \ / \ 

(2.1) Ipif)I,{g) = J2rl[P^j (^^j W2.(/®.5), 

as given, for instance, in Nualart's book [21], Proposition 1.1.3; the contrac- 
tion / 0r g is the element of '}{<^ip+<i-2r) defined by 

(/ 9){-^li ■ • ■ ) Sp—ri til • ■ • ) tq—r) 

(2.2) = / f{si,...,Sp-r,Ui,...,Ur) 

J[0,T]P+9-2'- 

X g{ti, . . . , tq-r, Ui, . . . , Ur) dui ■ ■ ■ dUr- 

We now introduce the Malliavin derivative for random variables in a chaos 
of finite order. If / G W®", we will use the following rule to differentiate in 
the Malliavin sense: 

DtlM)=nIn-iifn{;t)), tG[0,l]. 

It is possible to characterize the convergence in distribution of a sequence 
of multiple integrals to the standard normal law. We will use the following 
result (see Theorem 4 in [22], also [23]). 

Theorem 2.1. Fix n>2 and let {Fk,k > 1), F^ = In{fk) (with fk G 
"H®" for every k > 1), be a sequence of square-integrable random variables 
in the nth Wiener chaos such that E[F^] ^1 as k ^ oo. The following are 
then equivalent: 

(i) the sequence (i*A;)fc>o converges in distribution to the normal law 
Ar(0,l); 

(ii) B[F^] ^3 as k^oo; 
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(iii) for all 1 <l <n — 1 

(iv) IlL'Ffcll^ ^ n m ask 
tive with respect to B. 



it holds that limfe_,oo \\fk ®l /fc|l7^s}2{n-i) = 0; 

oo, where D is the Malliavin deriva- 



Criterion (iv) is due to [22]; we will refer to it as the Nualart-Ortiz-Latorre 
criterion. A multidimensional version of the above theorem has been proven 
in [24] (see also [22]). 

3. Variations for the Rosenblatt process. Our observed process is a Rosen- 
blatt process (-Z^(t))tG[o,i] with self-similarity parameter H G (i, 1). This cen- 
tered process is self-similar with stationary increments and lives in the sec- 
ond Wiener chaos. Its covariance is identical to that of fractional Brownian 
motion. Our goal is to estimate its self-similarity parameter H from discrete 
observations of its sample paths. As far as we know, this direction has seen 
little or no attention in the literature and the classical techniques (e.g., the 
ones from [5, 27] and [28]) do not work well for it. Therefore, the use of the 
Malliavin calculus and multiple stochastic integrals is of interest. 

The Rosenblatt process can be represented as follows (see [29] ) : for every 
tG [0,1], 

Z^{t):=Z{t) 



(3.1) 



diH) 



JO 



u,yi) 



X diK^' {u,y2)du 



dWiyi)dW{y2), 



where {W(t),t G [0, 1]) is some standard Brownian motion, is the stan- 
dard kernel of fractional Brownian motion of index H' (see any reference on 
fBm, such as [21], Chapter 5) and 

(2(2// - l))i/2 



(3.2) 



H' = ^^^ and d{H) 



For every t G [0,1], we will denote the kernel of the Rosenblatt process 



with respect to W by 

Lt{yi,y2) ■■=Lt{yi,y2) 

ft 



(3.3) 



■.= d{H) 



diK'''{u,yi)diK'''{u,y2)du 



l[0,i]2(?/l>2/2)- 



In other words, in particular, for every t, 

Z{t)=l2iLt{.)), 

where I2 denotes the multiple integral of order 2 introduced in Section 2. 
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Now, consider the filter a = { — 1, 1} and the 2-variations given by 

V (9 ^ 1 V (Zji/N) - Z{{z - 1)/N)y 
VN[Z,a) ^Z.E(Z(i/iV)-Z((z-l)/iV))2 



N 

2H 



1=1 



-n— 11 



The product formula for multiple Wiener-Ito integrals (2.1) yields 

hU? = W ®f) + 4l2(/ 01 /) + 2||/||i2([o,ip). 
Setting, for z = 1, . . . , N ^ 

(3.4) Ai := Li/N - L(j_i)/Ar, 

we can thus write 



^livj "^lirjj =ih{A^)f = hiA^A,)+4l2iA<S,iA,) + N~^'' 

and this implies that the 2-variation is decomposed into a fourth chaos term 
and a second chaos term: 

TV 

Vn{2, a) = N^^-^ Y.(^4iAi <S> Ai) + AhiA, ®i Ai)) := T4 + Tg. 
1=1 

A detailed study of the two terms above will shed light on some interesting 
facts: if H < J, then the term T4 continues to exhibit "normal" behavior 
(when renormalized, it converges in law to a Gaussian distribution), while 
the term T2, which turns out to be dominant, never converges to a Gaussian 
law. One can say that the second Wiener chaos portion is "ill behaved"; 
however, once it is subtracted, one obtains a sequence converging to AA(0, 1) 
for E (|, |), which has an impact on statistical applications. 

3.1. Expectation evaluations. 

3.1.1. The term T2. Let us evaluate the mean square of the second term, 

TV 

r2:=4iV2^^-i^/2(^i®i^i). 

4 = 1 

We use the notation /j = (^^, jj'] for i = 1, . . . , N . The contraction Ai Ai 
is given by 

{Ai^iAi){yi,y2) 
1 

Ai{x,yi)Ai{x,y2)dx 
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: d{H) 1^ dx 1[0A/N] {yi V x)l[oi/Ar] (1/2 V x) 
fi/N 

diK^' {u,x) diK^' {u,yi) du 

(3.5) - l[o,(i-i)/Af](yi Vx) 

r-(i-l)/7V 

diK^ {u,x)diK^ iu,yi)du 
diK" {v,x)diK" {v,y2)dv 

xVj/2 

I[0,(i-1)/Ar](y2 Vx) 

r-(i-l)/Af , , 

diK^ {v,x)diK" {v,y2)dv 
Defining 

(3.6) a{H) := H'{2H' - 1) = H{H + 1) /2, 
note the following fact (see [21], Chapter 5): 

(3.7) / diK'' {u,yi)diK'' {v,yi)dyi = a{H)\u-v\^''-^; 
Jo 

in fact, this relation can easily be derived from Jq^^ K^' {u,yi)K^' {v,yi) dyi - 
R^' (n, v) and will be used repeatedly in the sequel. 

To use this relation, we first expand the product in the expression for the 
contraction in (3.5), taking care to keep track of the indicator functions. The 
resulting initial expression for {Ai ®i ^i)(2/i;y2) contains four terms, which 
are all of the following form: 

Ca,b ■■= d{Hf dx l[o,a] {yi V x)l[o,6] (2/2 V x) 

X r diK^'{u,x)diK^'{u,yi)du 

J M=3/i Vx 
b 

diK^\v,x)diK^\v,y2)dv. 

v=y2\/x 

Here, to perform a Pubini argument by bringing the integral over x inside, 
we first note that x < u Av while u G [yi, a] and v £ [y2,b]. Also, note that 
the conditions x < u and u < a imply that x < a and thus l[o,a](2/i V x) 
can be replaced, after Fubini, by l[o,a](yi)- Therefore, using (3.7), the above 
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expression equals 

Ca,b = d{H)'^l[Q,,-^y,[Qf,-^{yi,y2) diK^ {u,yi)du diK^ {v,y2)dv 

Jyi •''1/2 



X 



ruAv ^ ^ 

/ diK" {u,x)diK^ {v,x)dx 
Jo 



= d{H)^aiH)l[o,a]xmiyi^y2) r f diK^'iu,y,)diK'''iv,y2) 

Ju=yi Jv=y2 

X |n — dudv 
= d{Hfa{H) r f diK{u,yi)diK^'{v,y2)\u-v\'^^'~^dudv. 

J u=yi J v=y2 

The last equality comes from the fact that the indicator functions in yi , y2 
are redundant: they can be pulled back into the integral over dudv and, 
therein, the functions diK^ {u,yi) and diK^ {v,y2) are, by definition as 
functions of yi and y2, supported by smaller intervals than [0,a] and [0,6], 
namely [0, u] and [0, v] , respectively. 

Now, the contraction {Aii^iAi){yi,y2) equals Cj/jy^j/Ar + C(j_i)/Ar^(i_i)/Ar - 
C{i-i)/N,i/N ~ Q/Ar,{i-i)/Af • Therefore, from the last expression above, 

{Ai(^iAi){yi,y2) 

/ pi/N ri/N 

= a{H)d{Hf[ du dvdiK^ {u,yi) 
\Jyi Jy2 

xd,K^'{v,y2)\u-v\'^'-' 

ri/N r{i-l)/N 



/ du dvdiK" {u,yi) 

Jyi •'y2 

H' f„. „, M„, „.\2H'-2 



(3.8) xdiK"" {v,y2)\u-v 

i-(i-l)/N fi/N 



/ du dvdiK" {u,yi) 

Jyi Jy2 

xd^K'''{v,y2)\u-v\^'''-^ 

r-{j-l)/Af /•{i-l)/N 



+ / du dvdiK^ {u,yi) 

Jyi J s/2 

xd,K^\v,y2)\u-v\'''''-^y 



Since the integrands in the above four integrals are identical, we can simplify 
the above formula, grouping the first two terms, for instance, to obtain an 
integral of v over Jj = (^^, jj]-, with integration over u in [yi, The same 
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operation on the last two terms gives the negative of the same integral over 
V, with integration over u in [yi,^^]. Then, grouping these two resulting 
terms yields a single term, which is an integral for {u,v) over li x li. We 
obtain the following, final, expression for our contraction: 

(3.9) 

= a{H)d{Hf [ [ diK^' {u,yi)diK"' {v,y2)\u-v\'^^'-^ dudv. 

J JliXli 

Now, since the integrands in the double Wiener integrals defining T2 are 
symmetric, we get 

N 

E[r|] = N^^~h6 ■ 2! J2 ®i ®i ^i)L2([o,i]2). 

To evaluate the inner product of the two contractions, we first use Fubini 
with expression (3.9); by doing so, one must realize that the support of 
diK^ {u,yi) is {u > yi}, which then makes the upper limit 1 for the inte- 
gration in yi redundant; similar remarks hold with respect to u',v,v' and 
2/2- In other words, we have 

{Ai (g)i Ai,Aj ^j)i2([o,i])2 

= a{Hfd{H)'^ [\y1dy2 [ III du'dv'dudv 
Jo Jo Jli Jh Jlj Jlj 

X \u-v\^^'-^\u -v\^^'-^diK^\u,yi)diK^\v,y2) 

xdiK'''{u',yi)diK'''{v',y2) 

(3.10) =a{Hfd{Hf[ j j j \u-v\'''''-''\v! -v'\^'''-''du'dv'dvdu 

X / diK"" {u,yi)diK'' {u',yi)dyi 
Jo 

rvAv' 

X / diK"" {v,y2)diK'' {v',y2)dy2 
Jo 

= a{Hfd{Hf( III |n-i;|2^'-V-T;f^'->-nf^'-2 
Ju Jh Jlj Jlj 

X \v — V \ ~ du dv dvdu, 

where we have used the expression (3.7) in the last step. Therefore, we 
immediately have 

E[r|] = N^^-'^2>2a{Hfd{Hf 
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N 

(3.11) X y] / / / / du'dv'dvdu 

^, I \2H'-2\ I /\2H'-2\ /\2H'-2\ /\2H'-2 

X \u — v\ \U — V \ \u — u\ \V — V \ 

By Lemma 5 in the Appendix, we conclude that 

Jim EiniA^-- = 6MHyw{^^ - ^ 

(3.12) = l6d{Hf 



3.1.2. The term T4. Now, for the L^-norm of the term denoted by 

N 

^4:=7V2^^-^^/4(^.®^^), 
1=1 

by the isometry formula for multiple stochastic integrals, and using a correc- 
tion term to account for the fact that the integrand in T4 is nonsymmetric, 
we have 

N 

E[r|] = 8iV4^-2 J2 {Ai ® A; A, ® Aj) l2([o,i]^) 

N 

+ 4A^4H-2 J2 4(Ai®i^j;^i®i^i)L2{[o,i]2) =-%,o + %,i- 

We separate the calculation of the two terms 74^0 and 74^1 above. We will see 
that these two terms are exactly of the same magnitude, so both calculations 
must be performed precisely. 

The first term, 74,0) can be written as 

N 

'?4,o = 8iV^^~^ ^ |(^j,ylj-)^2([o,i]2)| . 
We calculate each individual scalar product (^j, Aj)j;^2([o,i]2) as 

(^i,^j)L2{[0,l]2) 
»1 /-l 

Ai{yi,y2)Aj [yi , ^2) dyi dy2 



Jo 



\2 



d{H) dyidy2l[o,i/N/\j/N]{yi\/ y2) 
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X / (u,yi)9iK^ (u,y2)d^i- l[o,{i-i)/JV](yi Vy2) 



(i-l)/JV 



diK''\u,yi)diK''\u,y2)du 



X / dxK^ {v,yi)diK" (^^, ^2) - l[o,0-i)/iV] (yi V ^2) 
\Jy1Vy2 



0-l)/Af 



diK''\v,yi)diK'''{v,y2)dv 



H', 



i/N j-j/N 



d{Hf 

l{i-l)/N J{j-1)/N 

Here, (3.7) yields 



yiVj/2 

uAv 



du dv 



diK'''{u,yi)diK"'iv,yi)dy 







{A,A,)L2^[0,m=d{Hfa{H)^ f I 



\u — v\'^^ "^dudv. 



where, we have agam used the notation Ij = {^j^, jf] for i = 1, . . . ,N . We 
finally obtain 



(3.13) 



d{Hfa{Hf 1 
H{2H-1) 2 



2 


i-j 


2H 




2H 




2Hi 




N 




N 




N 





where, more precisely, d{Hfa{Hf{H{2H - 1))"^ = 2. Specifically, with the 
constants ci^h-, C2,h and c'^ jj given by 



ci 



,j^:=2 + £(2A;2^^-(fc-l)2^-(A: + l 



,2H\2 



k=l 



(3.14) 



C2,H:=2F2(2i7-l)/(4//-3), 

c',^H-=i'^H{2H-l)f = 9/16, 

using Lemmas 3, 4 and an analogous result for H = 3/4, we get, asymptot- 
ically for large N, 

(3.15) jim A^r4,o = 16ci,iy, l/2<if<|, 



(3.16) 
(3.17) 



W— >oo 

4-4_ff/7- _ -ic„ u ^ 3 



4' 



lim iV*-*^r4o = 16c2H, H> 
lim —74 = 16c'i rr = 16, H = 



The second term, 74^1, can be dealt with by obtaining an expression for 
{Ai (g)i Aj;Aj 01 ^i)i2([o_i]2) 
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in the same way as the expression obtained in (3.10). We get 

N 

r4,i = 16N^^-^ J2 i^i ®i ®i ^i) LH[0,1]^) 

N rl rl rl rl 

I ,1 I 



^0 



dydzdy dz 



x\y- z + i- j\^^'~'^\y' - z' + i- ip^'"^ 
x\y-y' + i- jf -z' + i- jf 

Now, similarly to the proof of Lemma 5, we find the the following three 
asymptotic behaviors: 

• if G |), then TijjNT4^^i converges to 1, where 

(3.18) ri,H := 16d{H)^a{Hyci,H; 

• if H > J, then r^^A^^~^^74^i converges to 1, where 

(3.19) T2H-=^2d{Hfa{Hf f\l - x)x'^"-^ dx; 

Jo 

• if H = J, then r^^(iV/ log iV)74^i converges to 1, where 

(3.20) T3,// := 2,2d{Hfa{Hf. 

Combining these results for 74^1 with those for 74^0 in lines (3.15), (3.16) 
and (3.17), we obtain the asymptotics of E[T4] as — > co: 

hm iVE[T|] = ei,H, ifi?G(i,|); 

A' — ^00 

hm iV4-4HE[r2] ^ if G (1, 1); 

N 3 

where, with u^h, ^ = 1,2,3, given in (3.18), (3.19) and (3.20), we defined 

(l/2)ci,H+ri,H, 

il/2)c2,H+T2,H, 
C3,H + r3 H- 



(3.21) e2,H 

e3,H 



Taking into account the estimations (3.15), (3.16) and (3.17), with c^^h 
in (3.12), we see that E[T4] is always of smaller order than E[r2]; therefore. 
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the mean-square behavior of V/v is given by that of the term T2 only, which 
means that we obtain, for every H > 1/2, 



(3.22) hm E 



3.2. Normality of the fourth chaos term T4 when H < 3/4. The calcula- 
tions for r4 above prove that limAr^oo E[G^] = 1 for < 3/4, where ei^n is 
given in (3.21) and 

(3.23) Gn ■■= ViViV2H-ig-i/2^^ ^ ^.^ _ 

Similarly, for i7 = |, we showed that limAr_+oo ^[(5^] = 1, where e^^n is given 
in (3.21) and 



(3.24) := ^^N'H-'e^},h ® • 

Using the criterion of Nualart and Ortiz-Latorre [part (iv) in Theorem 2.1], 
we prove the following asymptotic normality for Gn and G^. 

Theorem 3.1. If H e (1/2,3/4), then Gn given by (3.23) converges in 
distribution as 

(3.25) lim GN=N{<d,l). 

AT— »oo 

If H = 3/4, then Gn given by (3.24) converges in distribution as 

(3.26) lim Gn =M{0,1). 

Proof. We will denote by c a generic positive constant not depending 
on N. 

Step (Setup and expectation evaluation). Using the derivation rule 
for multiple stochastic integrals, the Malliavin derivative of Gn is 

TV 

DrGN = VNN^^'-^e-^ihY^hiiAi ® Ai){-,r)) 

i=l 

and its norm is 



l-DG'Af|li2([o,i]) 



N ,1 



iV^^-Heer}^ E / drh{{A®Ai){;r))h{{Aj^Aj){;r)). 
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The product formula (2.1) gives 

II^G'Af|li2([0,l]) 

= iV4^-il6eri 



N 1 

X V / dr[lQ{{Ai0Ai){-,r)^{Aj(^Aj){-,r)) 
■ 1 Jo 



+ 9h{iAi Ai){-,r) 01 (Aj (g) Aj){-, r)) 

+ 9l2{{Ai (g Ai){-,r) 02 {Aj Aj){; r)) 

+ 3\Io{{Ai^Ai)i-,r) 03 {Aj^Aj){-,r))] 

=: Jg + J4 + J2 + </o- 

First, note that, for the nonrandom term Jq that gives the expected value 
of the above, we have 

N 

Jo = 16er]^A^^^-'3! ^ / 2/2)^^(^3, 2/4)^,(2/1, ^2) 

X Aj (2/3 , 2/4) ^2/1 c?2/2 ^2/3 c?2/4 

N 

This sum has already been treated: we know from (3.15) that Jo/4 con- 
verges to 1, that is, that lim^^^ E[||L)GAr|||2([o 1])] = 4. This means, by the 
Nualart-Ortiz-Latorre criterion, that we only need to show that all other 
terms Je, J4, J2 converge to zero in L'^{Q) as N ^ 00. 

Step 1 (Order-6 chaos term). We first consider the term Jg: 

N 1 

J6 = ciV4^-iV / drlG{{Ai®Ai){;r)^{Aj^Aj{;r))) 

N 

= cN^"-^ J2 ^6((^* ® A,) {A, 01 Aj)). 

We study the mean square of this term. We have, since the L^-norm of 
the symmetrization is less than the L^-norm of the corresponding unsym- 
metrized function, 

B ^eiiAi Aj) {Ai 01 Aj))] 



VARIATIONS AND ESTIMATORS VIA MALLIAVIN CALCULUS 



19 



< 6! ^ {{Ai (g) Aj) (g) {Ai 01 Aj), {Ak Ai) ® {Ak ®i Ai)) Lmo,i]^) 

i,j,k,l 

= 6! ^ (Aj, Afc)i2([o^l]2)(Aj, Az)^2([o,l]2)(^j ®l Aj,Ak 01 Ai)i2([o^l]2). 
i,j,k,l 



We get 



E[J|] < cA^«^"2 y [ du( dv I du' I dv' 
ij^k/i. Jij Jh Jh 



X 


2 


i — k 


2H 


i-k + 1 


2H 


i - 


-k-l 


2H 


N 




N 






N 




X 


2 


j-l 


2H 


j-l + l 


2H 


3 - 


-l-l 


2H-I 




N 




N 






N 




First, we show that for H G (1/2,3/4), 


we have. 


for large 


N, 





(3.27) 



E[J|] < cN' 



With the notation as in Step 1 of this proof, making the change of variables 

irly for the ( 

1 1 1 



u = {u — ^^)N, and similarly for the other integrands, we obtain 



E[J|] < cN 



m~2 



i,j,k,l 



du dv du dv' 



[0,1]* 

X \u-v + i-j\^"'~'^\u-u' + i-k\'^^'-'^ 
^\u'-v + j- A:|2^'-> -v' + k- 
x(2|i-fc|2^-K-/fc + lp^-K-fc-lp^ 

xi2\j-ir-\j-i+ir-\j-i-in 



-y 



i,j,k,l 



du dv du dv' 



x\u-v + i- V -u + i- 

X |n' - + j - fc^^'-^l^; -v' + k- Zp^'^^ 

X {2\i - -\i-k + 1|2^ -\i-k- 1|2^) 

x(2|j-Zr-|,-Z + l|2^-U-Z-l|2^). 

Again, we use the fact that the dominant part in the above expression is 
the one in where all indices are distant by at least two units. In this case. 



20 



C. A. TUDOR AND F. G. VIENS 



up to a constant, we have the upper bound \i — k]^^^"^ for the quantity 
{2\i — fcp^ — \i — k + l]"^^ — \i — k — Ip^)- By using Riemann sums, we can 
write 



1 



1 



i,j,k,l 



i j k I 



^8H'-8jy4-H'-4 



where / is a Riemann integrable function on [0, 1]^ and the Riemann sum 
converges to the finite integral of / therein. Estimate (3.27) follows. 



Step 2 (Chaos terms of orders 4 and 2). To treat the term 

N 1 

J4 = cA^4H-i^ / drh{{Ai(^Ai){;r)(^i{Aj(^Aj){;r)), 

since 14(g) =I'i{g), where g denotes the symmetrization of the function g, 
we can write 

N 

J4 = cN^^"^ J2 {Ai,Aj)L2^Q^i]2h{Ai^Aj) 



N 



+ cN"^^ {Ai ®i Aj) ® {Ai ®i Aj) 

=: Ji^i + J4^2- 

Both terms above have been treated in previous computations. To illustrate 
it, the first summand J^^i can be bounded above as follows: 

N 

E| J4,l|^ < CA^®-^"^ Y (Aj, Aj)^2([o,l]2)(Aj,^fc)L2([0,l]2) 
i,j,k,l=l 

X (^fc> A)L2([0,l]2)(^j, A)l2([o^1]2) 



N 

i,j,k,l=l 



i-j + l 
N 



2H 



+ 



N 



2H 



N 



2H- 



+ 



i-k + l\^" 

N ) 



N 

N 



N ) 
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+ 



N J \ N 

and, using the same bound c\i — Jp^~^ for the quantity |i — j + + — 
j — Ip-'^ — 2|i — jl"^^ when \i — j\ > 2, we obtain 

N 

i,j,k,l=l 

^^^SH^e± f \^-J\'''-'\^-kr-'\J-l?''-'\k-l\'^-' 

i,j,k,l=l 

This tends to zero at the speed A^^H-e as ^ oo by a Riemann sum argu- 
ment since H < j. 

One can also show that E|J4^2p converges to zero at the same speed 
because 

TV 

i,j,k,l=l 

N 



{ {\u-v + i-j\ 



x\u -v' + k-l\ 
X \u — u + i — k\ 



X \v-v' +j-l\f"'-'^dv'dudvdu^ 



Thus, we obtain 

(3.28) E[J|] < cN^^-^ 

A similar behavior can be obtained for the last term J2 by repeating the 
above arguments: 

(3.29) E[J|] < cN^^-^. 

Step 3 (Conclusion). Combining (3.27), (3.28) and (3.29) and recalling 
the convergence result for E[T|] proven in the previous subsection, we can 
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apply the Nualart-Ortiz-Latorre criterion and use the same method as in 

3 
4 



the case < | for = 3/4, to conclude the proof. □ 



3.3. Nonnormality of the second chaos term T2 and limit of the 2-variation. 
This paragraph studies the asymptotic behavior of the term denoted by T2 
which appears in the decomposition of Vj\f{2, a). Recall that this is the dom- 
inant term, given by 



and, with ^c^^h = ^d{H) given in (3.12), we have shown that 



^lim^E[(iVi-^r2C3-^^)2] = l. 



With Tat := N^''^T2C^'h \ one can show that in L\Q), 

Jim^||Dr^||i2([o,i])=2 + c, 

where c is a strictly positive constant. As a consequence, the Nualart-Ortiz- 
Latorre criterion can be used to deduce that the T/v do not converge to the 
standard normal law. However, it is straightforward to find the limit of T2, 
and thus of VXr, in L^(r2), in this case. We have the following result. 

Theorem 3.2. For all H £ (1/2,1), the normalized 2-variation N^~^VNi2, 
a)/{Ad{H)) converges inL'^{Q) to the Rosenblatt random variable Z{1). Note 
that this is the actual observed value of the Rosenblatt process at time 1. 

Proof. Since we already proven that N ~ converges to in L2(n), 
it is sufficient to prove that N'^-^T2/{Ad{H)) - Z{1) converges to in ^^(O). 
Since T2 is a second-chaos random variable, that is, is of the form l2ifN), 
where /at is a symmetric function in L^([0, l]'^), it is sufficient to prove that 

converges to Li in L2([0, 1]2), where Li is given by (3.3). From (3.9), we get 
fN{yi,y2) = ^N^''~'aiH)d{Hf 



(3.30) X J2 



N 



i=l 



X diK^' {u,yi)diK"' {v,y2)dudv]. 
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We now show that jj(^fN converges pointwise, for yi,y2 G [0, 1], to the 
kernel of the Rosenblatt random variable. On the interval li x li, we may 
replace the evaluation of diK^ and diK^ at u and v by setting u = v = 
i/N. We then get that /Ar(yi,y2) is asymptotically equivalent to 

N 

AN^''-\{H)d{Hf ^ h,N>y,yy, diK^"' {i/N, y{) d^K^'' {i/N, y^) 
1=1 

X / / dudv\u-v\^^'~'^ 

J JliXli 

1 ^ 

= AN''-'d{Hf-Y.^i/N>yiyy29iK'''{i/N,yi)diK'''{i/N,y2), 

i=l 

where we have used the identity JJj^xii dudv \u — v\'^^'~'^ = a{H)~^ N~'^^' = 
a{H)-^N-^-^. Therefore, we can write, for every yi,?/2 G (0, 1)^, by invok- 
ing a Riemann sum approximation, 

Af->oo M{H) 

1 ^ 

= d{H) \\-m^ - h/N>y,yy, {i/N, yi) diK'''{i/N, ya) 

i=l 

= d{H) f' 9iK^'(n,yi)aii^^'(u,y2)dn = Li(yi,y2). 

To complete the proof, it suffices to check that the sequence (4(i(if))~^ x 
N^-^/n is Cauchy in ^^([0, l]^) [indeed, this implies that {M{H))-^ N^'^ Jn 
has a limit in -L^([0, 1]^), which obviously coincides with the a.e. limit Li and 
then the multiple integral l2{{4:d{H))~^N^~^ f]\f) will converge to l2{Li)]. 
This can be checked by means of a straightforward calculation. Indeed, one 
has, with C{H) a positive constant not depending on M and A^, 

In - ^^~'^/M|li2([o,i]2) 



N 

C{H)N'^ I 11 I 
ij^iJiiJh-JijJij 



I 12^-21 / /\2H'-2 

\u — v\ \U — V 



X In - u'l^'^'-^lv - v'l^'^'-^du'dv'dudv 



+ C{H)M'^^ 

^ njM fi/M fj/M rj/M 



^J(i-l)/MJ{i-l)/M J{j-l)/M J{j-1)/M 



x^l' I' 1^' 1^' |n-T;|2^'-2|^/_^/,2/f'-2 
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t\2H'-2 



x\u — u 



x\v-v'\'^^ dv! dv' dudv 



l-H tktI-H t^j2H-1^2H-1 



2C{H)M'-''N'-"M 

N M j/M 

-EE 

i=ij=i 



IiJliJ(j~l)IM J{j~l)/M 



j/M 



du' dv' du dv 



X |n-n'P^'-2h-7;'|2^^'-2. 



The first two terms have already been studied in Lemma 5. We have shown 
that 



N 

n'^'j: f f f f 



X In - n?^ - ^;'|2^^ -2 du' dv' dudv 



converges to {a{H)'^ H {2H — \))~^ . Thus, each of the first two terms in (3.31) 
converge to C{H) times that same constant as M,N go to infinity. By the 
change of variables which has already been used several times, u= {u — , 
the last term in (3.31) is equal to 



C{H){MN) 



H 



1 



{NM) 



2H'~2 



N M 

EE 



du dv du dv' 



C{H) 
MN 



\J[0,1]^ 








X 


u — 


vl^"" 




X 


U 

N 


u' 
M 




2H'- 


X 


V 


v' 




2H'~ 


N 


M 




N M 

EE 

i=lj=l 


1 

J[0,l]4 


du dv du' dv' 



^\^_^\2H'-2\J _n2H'-2 
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V v' i j 
^ N~M^N~M 

For large the term — ^ in front of — is negligible and can be ig- 
nored. Therefore, the last term in (3.31) is equivalent to a Riemann sum than 
tends, as M,N oo, to the constant (/q /q \u — v\'^^'~'^ dudv)"^ /q \x — 
y|2{2H'-2)_ rpj_^jg precisely equal to 2{a{HfH{2H-l))-^, that is, the limit 
of the sum of the first two terms in (3.31). Since the last term has a leading 
negative sign, the announced Cauchy convergence is established, completing 
the proof of the theorem. □ 



Remark 1. One can show that the 2- variations VAr(2,a) converge to 
zero almost surely as N goes to infinity. Indeed, the results in this section 
already show that VAr(2, a) converges to in L^(0), and thus in probability, 
as N ^ oo; the almost sure convergence is obtained by using an argument 
in [4] (proof of Proposition 1) based on Theorem 6.2 in [6] which gives the 
equivalence between the almost sure convergence and the mean-square con- 
vergence for empirical means of discrete stationary processes. This almost- 
sure convergence can also be proven by hand in the following standard way. 
Since V]s[{2,a) is in the fourth Wiener chaos, it is known that its 2pth mo- 
ment is bounded above by Cp(E[{Vn{2, a))^])^/^, where Cp depends only on p. 
By choosing p large enough, via Chebyshev's inequality, the Borel-Cantelli 
lemma yields the desired conclusion. 

3.4. Normality of the adjusted variations. According to Theorem 3.2, 
which we just proved, in the Rosenblatt case, the standardization of the 
random variable Vat (2, a) does not converge to the normal law. But, this 
statistic, which can be written as Vn = T4 + T2, has a small normal part, 
which is given by the asymptotics of the term T4 , as we can see from Theorem 
3.1. Therefore, Vjv — T2 will converge (under suitable scaling) to the Gaussian 
distribution. Of course, the term T2, which is an iterated stochastic integral, 
is not practical because it cannot be observed. But, replacing it with its limit 
Z{1) (this is observed), one can define an adjusted version of the statistic 
Vn that converges, after standardization, to the standard normal law. 

The proof of this fact is somewhat delicate. If we are to subtract a mul- 
tiple of Z{1) from Vat in order to recuperate T4 and hope for a normal 
convergence, the first calculation would have to be as follows: 

Vn{2, a) - ^Z{1) = VV(2, a)-T2 + T2- ^Z{1) 



(3.32) 



T 



-H 



T2 - Z{1) 
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The term Ta, when normahzed as Ta, converges to the standard normal 
law, by Theorem 3.1. To get a normal convergence for the entire expression 
in (3.32), one may hope that the additional term U2 ■= ^i^'h [^g|=g^2 — -^(1)] 
goes to "fast enough." It is certainly true that U2 does go to 0, as we just 
saw in Theorem 3.2. However, the proof of that theorem did not investigate 
the speed of this convergence of U2 - For this convergence to be "fast enough," 
one must multiply the expression by the rate ^/N which is needed to ensure 
the normal convergence of T4: we would need f/2 <C N~^/'^ . Unfortunately, 
this is not true. A more detailed calculation will show that U2 is precisely 
of order \/iV. This means that we should investigate whether ^/NU2 itself 
converges in distribution to a normal law. Unexpectedly, this turns out to 
be true if (and only if) H <2/3. 

Proposition 2. With U2 as defined in (3.32) and H < 2/3, we have 
that \/iVf^2 converges in distribution to a centered normal with variance 
equal to 

(3.33) ■.= 32d{HMHfp^k'''-'FQ-^, 

where the function F is defined by 

F{x)= I dudvdu dv'\{u-u)x + l\^^'-'^ 

X [a{H)Wu - v\\u' - v'\\{v - v')x + l\f"'-^ 

(3.34) 

- 2aiH){\u -v\\iv- u')x + l\f^'-^ 

+ \^u-u')x + l\'^'-% 

Before proving this proposition, let us take note of its consequence. 

Theorem 3.3. Let {Z{t),t G [0,1]) be a Rosenblatt process with self- 
similarity parameter H € (1/2,2/3) and let previous notation for constants 
prevail. Then, the following convergence occurs in distribution: 



lim 

]V-,cx, ^ei^H + fl,H 



V^{2,a)-^Z{1) 



:AA(0,1). 



Proof. By the considerations preceding the statement of Proposition 
2, and (3.32) in particular, we have that 

V^{2,a)-^Z{1] 



NTi + ^NU2 
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Theorem 3.1 proves that \'NT4 converges in distribution to a centered nor- 
mal with variance ei^n- Proposition 2 proves that y/NU2 converges in distri- 
bution to a centered normal with variance fi^H- Since these two sequences 
of random variables live in two distinct chaoses (fourth and second respec- 
tively), Theorem 1 in [24] implies that the sum of these two sequences con- 
verges in distribution to a centered normal with variance ei^n + fi,H- The 
theorem is proved. □ 

To prove Proposition 2, we must first perform the calculation which yields 
the constant fi^n therein. This result is postponed to the Appendix, as 
Lemma 6; it shows that E[(\/iVC/2)^] converges to fi^H- Another (very) 
technical result needed for the proof of Proposition 2, which is used to guar- 
antee that y/NU2 has a normal limiting distribution, is also included in the 
Appendix as Lemma 7. An explanation of why the conclusions of Propo- 
sition 2 and Theorem 3.3 cannot hold when H >2/3 is also given in the 
Appendix, after the proof of Lemma 7. We now prove the proposition. 

Proof of Proposition 2. Since U2 is a member of the second chaos, 
we introduce notation for its kernel. We write 



N 

— U2 = h{gN), 

\J Jl,H 

where ^at is the following symmetric function in L^([0, 1]^): 

9N{yi,y2) ■■= \ u{H) ^^^^^'^'^^ ~ -^I'^^i'^^) y 

Lemma 6 proves that E[(/2((7Ar))^] = Hfi^Af ||^2([o i]2) converges to 1 as ^ 00. 
By the result in [23] for second-chaos sequences (see Theorem 1, point (ii) 
in [23], which is included as part (iii) of Theorem 2.1 herein), we have that 
12(9 n) will converge to a standard normal if (and only if) 

jV™o llffiV ®1 S'Ar|li2([o,i]2) = 0, 

which would complete the proof of the proposition. This fact does hold if 
H <2/3. We have included this technical and delicate calculation as Lemma 
7 in the Appendix. Following the proof of this lemma is a discussion of why 
the above limit cannot be when H > 2/3. □ 

4. The estimators for the self-similarity parameter. In this section, we 
construct estimators for the self-similarity exponent of a Hermite process 
based on the discrete observations of the driving process at times 0, . . . , 1. 
It is known that the asymptotic behavior of the statistics V/v(2,o) is related 
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to the asymptotic properties of a class of estimators for the Hurst parameter 
H. This is mentioned in, for instance, [4]. 

We recall the setup for how this works. Suppose that the observed process 
X is a Hermite process; it may be Gaussian (fractional Brownian motion) or 
non-Gaussian (Rosenblatt process, or even a higher order Hermite process). 
With a = {—1, +1}, the 2-variation is denoted by 

(4.1) 5„(2,a) = l£(^(A)_x(l^ 

1 = 1 

Recall that E[S'Ar(2,a)] = N^"^^ . By estimating E[5Ar(2,a)] by 5Ar(2,a), we 
can construct the estimator 

iA o\ fo \ log Sn {2, a) 

(4.2) H^{2,a) = - ^^^^^ , 

which coincides with the definition in (1.4) given at the beginning of this 
paper. To prove that this is a strongly consistent estimator for H, we begin 
by writing 

l + VN{2,a)=SN{2,a)N^'', 

where Vn is the original quantity defined in (1.3), and thus 

log(l + Vn{2, a)) = log 5jv(2, a) + 2H\og N = -2{Hn{2, a) - H) log N. 

Moreover, by Remark 1, VN{2,a) converges almost surely to and thus 
log(l + VAr(2,a)) = VAr(2,a)(l + o(l)), where o(l) converges to almost 
surely as N ^ oo. Hence, we obtain 

(4.3) VN{2,a) = 2{H - Hm{2, a))(log Af)(l + o(l)). 

Relation (4.3) means that the Viv's behavior immediately gives the behavior 
of Hn - H. 

Specifically, we can now state our convergence results. In the Rosenblatt 
data case, the renormalized error — H does not converge to the normal 
law. But, from Theorem 3.3, we can obtain an adjusted version of this error 
that converges to the normal distribution. 

Theorem 4.1. Suppose that H > ^ and that the observed process Z is 
a Rosenblatt process with self-similarity parameter H. Then, strong consis- 
tency holds for Hn, that is, almost surely, 

(4.4) lim HN{2,a) = H. 

In addition, we have the following convergence in L^(i7); 

(4.5) hm ^l^log{N){HN{2,a)-H) = Z{l), 
Ar->oo 2d{H) 



VARIATIONS AND ESTIMATORS VIA MALLIAVIN CALCULUS 



29 



where Z{1) is the observed process at time 1. 

Moreover, if H < 2/3, then, in distribution as A^— >oo, with cs^h, ei,// 
and fi^H in (3.12), (3.21) and (3.33), 



N 



\/ei,H + fl,H 



-21og(iV)(i?^(2,a)-i/)-^Z(l) 



.AA(0,1). 



Proof. This follows from Theorems 3.3 and 3.2 and relation (4.3). □ 



APPENDIX 

Lemma 3. The series Efcli(2A;2^ - {k - l)^^ - {k + l)^^)^ is finite if 
and only if H € (1/2,3/4). 

Proof. Since 2k'^^ -{k- if" -{k + if^ = k'^"f{l), with f{x) := 2 - 
(1 — x)^^ — (1 + x)^^ being asymptotically equivalent to 2H{2H — l)x'^ for 
small X, the general term of the series is equivalent to {2H)^{2H — 1)^A;^^~^. 

□ 

Lemma 4. When H G (3/4,1), iV' E„=i,...,iV;|i-i|>2(2|y - 
|i^|2H_ |l^|2//)2 converges to H^{2H - 1)/{H - 3/4) asN^oo. 

Proof. This is left to the reader. The proof can be found in the ex- 
tended version of this paper, available at http : / / arxiv . org/ abs/0709 . 3896v2. 

□ 

Lemma 5. For all H >l/2, with I-i = (^, j^], i = l,...,N, 



N 

2H'-2 



\U — V\ 



hm N^'^y I I I I 

(A.l) 



X \v — V 



' ' ^ du dv dv du 



■ 2a{Hy 



2H-1 2H 



Proof. We again refer to the extended version of the paper, online at 
http://arxiv.org/abs/0709.3896v2, for this proof. □ 

Lemma 6. With fi^n given in (3.33) and U2 in (3.32), we have 

hm E[( ViVC/2)'] = fi,H. 

N~>oo 
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Proof. We have seen that ^/c^ = Ad{H). We have also defined 

-To - Zil] 



Let us simply compute the L^-norm of the term in brackets. Since this 
expression is a member of the second chaos and, more specifically, since 
T2 = hifN) and Z(l) = /2(Li), where Jn [given in (3.30)] and Li [given in 
(3.3)] are symmetric functions in L^([0, 1]^), it holds that 



E 



T2-Z{1) 



4d{H) 
Ad{HY 



Jn - Li 



L2([0,112) 



Af||L2([0,l]2) 



(/Af,-^^l)L2([0,l]2) + Il^l|li2([0,l]2)- 



M{H) 

The first term has already been computed. It gives 

J^2-2H 

■||/Af||L2([o,l]2) 



Ad{H)' 



N~^"a\H)d'{H) 

N 



xE 



du dv du dv' 



[0,1]^ 



X {\u — v\\u — v'\\u — u + i — j\\v — v' + i — j[f'^ ^. 

By using the expression for the kernel Li and Fubini's theorem, the scalar 
product of /at and Li gives 



4d{H) 



(/Ar,-^l)L2{[0,l]2) 



1 rl 



JO 



dyi dy2 



N 



l-H 



4d{H) 



N 



fN{yi,y2)Li{yi,y2) 



N^a{Hfd{HfY.I I 

i=i 

N 

N^a{Hfd{Hf [ f dudv f du 



dudv / du {\u — v\\u — u \\v — u 



t\\2H'-2 



{\u — t;| |u — u\ \v — u 



t\\2H'-2 
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= N~^^a{Hfd{Hf 

^ / (\u — v\\u — u' + i — j\\v — u' + i — jl)"^^ ""^ dudv du' . 

Finally, the last term ll-^^i |||2([q can be written in the following way: 
||Li||i2([o^y2) = d{Hfa{Hfj^ \u - n'|2(2^^'-2) dudu' 



'[0,1]2 
N 

\u-u'\^^^^'~'^dudu' 



= d{Hfa{HfN-^''y I \u-u' + i-j\''^^'''-^Uudu'. 

One can check that, when bringing these three contributions together, the 
"diagonal" terms corresponding to i= j vanish. Thus, we get 



where F is the function we introduced in (3.34). 

This function F is of class on the interval [0, 1]. It can be seen that 



F(0)= / dudvdu'dv' 



'[0,1] 

X {a{Hf{\u - v\\u' - v'lf^'-'^ - 2a{H)\u -v\ + l) 
a{Hf( I \u-v\^'''-A^ -2a{H) / \u - v\^''' du dv + I 



0. 



Similarly, one can also calculate the derivative F' and check that -F'(O) = 
0. Therefore, F{x) = o{x) as x — > 0. To investigate the sequence aN '■- 
N-^ Y.k=ii^ - ^ - we split it into two pieces: 



N~l 



a^ = N-^Y.iN-k-l)k^"~'F(\ 
k=i ^'^ 

N-l 



k=l k=l 

--: bN + CN- 
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Since ^at is the partial sum of a sequence of positive terms, one only needs 
to check that the series is finite. The relation F{l/k) <C 1/fc yields that it is 
finite if and only if 2H — 3 < —1, which is true. For the term cat, one notes 
that we may replace the factor A; + l by A; since, by the calculation undertaken 
for bi\f, N~^J2k=i^'^^~'^'^(l) converges to 0. Hence, asymptotically, we 
have 

N~l / -, s oo 

k=i ^'^^ k=i 

which thus converges to 0. We have proven that limaAr = limftjv = J2T=i • 
F{j:), which completes the proof of the lemma. □ 

Lemma 7. Defining 
we have limTv^oo Wqn ®i ^A^IIiaQo i]2) ~ ^ provided H < 2/3. 

— 1/2 

Proof. We omit the leading constant , which is irrelevant. Using 
the expression (3.30) for /at, we have 

gNiyi,y2) = N^''-'/'diH)a{H) 



N 

X 

i=l 



TV 

E/ / diK'''{u,y^)diK'''{v,y2)\u-v\^'''-^dvdu 
~^ Jh Jh 



- Li{yi,y2)- 

Here, and below, we will be omitting indicator functions of the type 
l[o,(i+i)/Af] (l/i) because, as stated earlier, these are implicitly contained in 
the support of diK^' . By decomposing the expression for Li from (3.3) 
over the same blocks li x li as for /at, we can now express the contraction 
qn '^1 dN as follows: 

(57V ®i 9N)iyi,y2) = N^^-\An - 2Bn + Cjv), 

where we have introduced three new quantities, 

AN:=N^"d{Hfa{Hf 



N 

X 



A" 

E / I dvdu I I dv' du 



X \\u-v\-\u' -v'\-\v-v'\f^'-^ 



dvdu I du \\u — v\ ■ \u — vW'^^ ^ 
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x5iK^'(n,yi)5iK^'K,y2), 

B^:=N"a{Hfd{HfJ2j J 

i=l 

x5iK^'(n,yi)aii^^'(n',y2) 

N 

= N^a{Hfd{Hf y I I dvdu I du'[\^ 
x5iK^'(n,yi)5ii^^'(n',y2) 



\u-v\ - W -v\?^ "2 



and 



CN = d{Hfa{H)f f dvdudiK"'{u,yi)diK"'iv,y2)\ 
Jo JO 

N 

= d{Hfa{H)Y^ / / dvdudiK"'{u,yi)diK^'{v,y2)\u-v\ 



hJlj 



The squared norm of the contraction can then be written as 

Wqn ®i S'Af|li2([o,i]2) 

= iV^^"2(pAr|||2([o,i]2) +4||5Ar|||2([o,i]2) 

+ II^Af Ili2([0,l]2) - 4(AAr, SAr)i2([o,i]2) 

+ 2(^Ar,C7v)i2([o^i]2) - A{Bn-,Cn)l'^{[0,1]2))- 

Using the definitions of A^r, and Cat, we may express ah six terms above 
exphcitly. All of the computations are based on the key relation (3.7). 
We obtain 

II^Ar||i2([o,i]2) 

= N^"a{Hfd{HYa{Hf 

N 

X ^ / j dvdu j / dv' du' 



i,j,k,l=l^^ 



J -^3 



X / dudv du dv 
Jik Jik -III Jh 

X \\u — v\ ■ \u! — v'\ ■ \v — v'\ ■\u — v\- \u — v\ 

X \v-v'\-\u-u\-\u' -u'W^"'-"^ 

N^"a{Hfd{Hf ^ ^ 
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> / du dv du dv' du dv dv! dv' 



X \\v — v' + i — j\-\v — v' + k — l\ 



II 11/ 'II- -II "/ -l\\2H'-2 

X m — t) • m — V \\u — v\\u' — V \\ 



and 



x\u-u + i-k\-\v! -u + j - /|]^'^'~^, 

ll-^A'lli2([o,i]2) 

= N'^^a{Hfd{Hf 



N 

X 

i,j,k,l=l ' 



N 

/ dvdu du / dudv du' 
, ,_, J/j J li J ij J If, J If, Jii 



n II' II- -I I-' -II -II' -'n2II'-2 

X [|m — v\ ■ \u — v\\u — v\ ■ \u — v\ ■ \u — u\ ■ \u — u W 
N^^a{Hfd{Hf 

X 7 / du dv du du dv du 

X — ■\u — V + i — j\\u — v\ ■\u — V + k — l\ 
x\u-u + i-k\ ■\u' -u' + j - 



l|CA^Ili2([o,i] 



N 

= N'^^a{Hfd{H)^ [ f dvdu [ f dv' du' 

i,j,kd=l -^^^ ■'^3 ■'^^ ■'^l 

X [|n — v\ ■ \u' — v'\ ■ \u — u'\ ■ \v — v'\\^^ ~^ 
= N'''a{Hfd{Hfj^^^^,^^ 

N 

X / dudvdu'dv' 

X [\u — v + i — j\ -{u — v' + k — l\ 
x\u-u' + i-k\ -Iv-v' 
The inner product terms can be also treated in the same manner. First, 

{An,Bn)l'2{[o,i]'2) 
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3-ff„/ U\7 7/ TT\4 



N^"a{H)'d{H) 

N 

X 



N 

2_. / dudv / du' dv' / dudv du 
-ri , J/, J/, J I, Ju Jii. Jh. J I, 



X \\u — v\ • In' — f'l • \v — v'\ -lu — v\ • \u' — v\ -lu — u\ • In' — u 



N7 J^7{2H'~2) 
N 

X > / du dv du' dv' du dv du' 



i,j,k,l=l 



X [|n — v\ • \u' — v'\ ■ \v — v' + i — j\ - lu — v\ 

X \u' — V + k — l\ ■ \u — u + i — k\ ■ \u' — u + j 

and 

= N'^^a{Hfd{Hf 



N 

X 

i,j,k,l=l ' 



N 

/ dudv / du dv' / du dv 
t1 -, J/, 7/, J/, J/, Jh. J I, 



p. 11/ /II /II II 11/ n2 f~f 

X [|n — u| • |n — u I • |n — n I • |n — n| • |n — n| • |n — v\ 
N^''a{HfdiH)'±^^^±^ 

N 

X / du dv du' dv' du dv 

X [|n — f I • |n' — n'l • |n — n' + i — j\ 

X \u — u + i — k\ ■ \u — V + k — l\ ■ u — v]'^^ 



Finally, 

= N"a{Hfd{Hf 



N 

X 

i,j,k,l=l ' 



N 

/ dudv du' / du dv 
t1 1 J I, hi J I, Jh. Jh 



X [|n — v\ ■ \u — n| • |u — v| • |u — n| • |n' — vW^^ ^ 
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= N^a{Hfd{Hf ^ ^ 



TVS j^5(2H'-2) 
N 

X / du dv du' du dv 

X \\u — v\ ■ \u' — V + i — j\ ■ \u — V + k — l\ 
x\u-u + i-k\-\u -v + j - 

We now summarize our computations. Note that the factors d{H)^ and 
'W Ni(2Hi-2) common to all terms. We also note that any terms corre- 
sponding to difference of indices smaller than 3 can be shown to tend col- 
lectively to 0, similarly for other "diagonal" terms in this study. The proof 
is omitted. We thus assume that the sums over the set D of indices i,j,k,l 
in {1, . . . , N} such that |i — j|, |A; — /|, |« — A;| and \j — l\ are all at least 2. 
Hence, we get 

\\gN ®i fl'Af|li2([o,i]2) 

(A.2) 

\i-j\.\k-l\.\i-k\.\j-l\-^^'-^ 



X E 

{i,j,k,l)eD 



X G 



1111 



— k — V i — k^ j — I 
where the function G is defined for {x,y,z,w) G [1/2, 1/2]^ by 
G{x,y,z,w) 

= a{H)^ / du dv du dv' du dv du dv' 

J[0,1]S 

X [\u - v\ ■ \u' - v'\ ■\u-v\- \u' - v'W^^'""^ 
X [\{v-v')x + l\ ■ \{v-v')y + l\ 

X \{u-u)z + l\-\(u' -u')w + lW^^'-'^ 

+ ^a{H)^ j dudvdu' dudvdu' 
J[0,l]e 

X \\u — v\ ■ \u — v\ ■\{u' — v)x -|- 1| • — v)y + 1| 

X |(u - u')z + 1| • |(m' - u')w + l|]2^'-2 

+ a{H)'^ / dudvdu dvdv' 

"'[0,1]'* 
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X [\{u-v)x + l\ ■ \{u' -v')y + l\ 
x\{u- u')z + l\-\{v- v')w + 

— 4a{Hy / du dv du' dv' du dv du' 

J[o,iV 

X [\u — v\ • \u' — v'\ • \u — v\ • \{v — v')x + 1| 

X \{u'-v)y + l\ ■ \{u-u')z + l\-\{u' -u')w + l\]'^^'-^ 

+ 2a{Hf / dudvdu' dv' dudv 

J [0,1] 6 

X [\u — v\ • \u' — v'\ • \{v — v')x + 1\ -{{u — v)y + 1| 
X |(n - u')z + 1| • |(u' - v)w + 



Aa{Hf 



du dv du' du dv 



[0,1]5 



x\\u — v\-\{v — u')x + 1| • |(u — v)y + 1| 

x\{u-u)z + l\ - \{u' - v)w + 

It is elementary to check that G and all its partial derivatives are bounded 
on [—1/2, 1/2]^. More specifically, by using the identity 



Jo Jo 



In — u| dudv, 



we obtain 

G(0, 0, 0, 0) = aiH)^ + Aa{H)^ + a(if)^ - 4a(i/)^ + 2a{H)^ - Aa{Hf 
= 0. 

The boundedness of G's partial derivatives implies, by the mean value the- 
orem, that there exists a constant K such that, for all {i,j,k,l) G D, 



G 



1111 



i — j k — I i — k j — I 



K K K K 

< ^ ^ + 7-, TT + T-. 77 + 



\k-l\ \i-k\ 

Hence, from (A. 2), because of the symmetry of the sum with respect to the 
indices, it is sufficient to show that the following converges to 0: 



(A.3) S:=N 



1 



E 



{i,j,k,l)&D 



\i- j\ -{k-ll-li- k\- \ j 

Ar4 



1 

K -il 
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We will express this quantity by singling out the term i' ■.= i — j and summing 
over it last: 

_ r,;\r3H-2 \ " r-i\H~2 ^ f \k — l\ ■ l^' + j — ^ ' \ j - ^ 

i'=3 (^i'+jj,k,l)&D ^ ^ 

For fixed i\ we can compare the sum over j, A;, / to a Riemann integral since 
the power H — 1 > —1. This cannot be done, however, for rather, 
one must use the fact that this is the term of a summable series. We get 
that, asymptotically for large N, 

N-l 
i'=3 

where the function g is defined on [0, 1] by 

rl — x pi rl 

(A. 4) g{x) := / / / dydzdw\(z-w){x + y-z){y-w)\^~^. 
Jo Jo Jo 

It is easy to check that g is a bounded function on [0,1]; thus, we have 
proven that for some constant K > 0, 

oo 
i'=3 

which converges to provided H < 2/3. This completes the proof of the 
lemma. □ 

We conclude this appendix with a discussion of why the threshold H <2/3 
cannot be improved upon, and the consequences of this. We can perform a 
finer analysis of the function G in the proof above. The first and second 
derivatives of G at = (0, 0, 0, 0) can be calculated by hand. The calculation 
is identical for dG/dx{0) and all other first derivatives, yielding [via the 
expression used above for a{H)], 

1 dG,-, 

= a{H)^ / dudv du dv' {v — v')[\u — v\ ■ \u — v'\]^~^ 

J[0,l]'i 

+ 4a{Hf [ dudvdu{v-u)\u-v\^-^ 
Jlo^r 
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+ a{Hf [ dudv{u-v) 



[0,1] 



-Aa{Hf / dudvdu' dv'{v-v')[\u-v\-\u' -v'W^~^ 
+ 2a{Hf I dudvdu' dv'[v-v')[\u-v\-\u' -v'\f-^ 

-Aa{Hf I dudvdu\v-u')\u-v\^-^. 

We note that the two hnes with Aa{H)^ cancel each other out. For each of 
the other four Unes, we see that the factor [v — v') is an odd term and the 
other factor is symmetric with respect to v and v' . Therefore, each of the 
other four factors is zero individually. This proves that the gradient of G at 
is null. Let us find expressions for the second derivatives. Similarly to the 
above calculation, we can write 



{l-H){2-H) 9x2 



■ a{Hf I dudvdu'dv'{v-v'f[\u-v\-\u'-v'\f-^ 

^[0,l]4 

+ 4a{Hf [ dudvdu{v-uf\u-vf-^ 
J[o,ip 



+ a(i/)^ / dudv{u-v)' 



'[0, 

-4aiHf f dudvdu'dv'{v-v'f[\u-v\-\u'-v'\f-^ 
+ 2a{Hf I dudvdu'dv'{v-v'f[\u-v\-\u'-v'\]"-^ 

J[0,l]4 

-Aa{Hf I dudvdu'iy -u'f\u-vf-^. 
i[o,ip 

Again, the terms with a{H)^ cancel each other out. The three terms with 
a{H)^ add to a nonzero value and we thus get 

{l-H){2-H) dx^ ^ ^ 

= -a{Hf I dudvdu' dv'{v-v'f[\u-v\-\u' -v'W^-^ 

J[0,l]4 



+ a{Hf [ dudv{u-vf. 
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While the evaluation of this integral is nontrivial, we can show that for all 
H > 1/2, it is a strictly positive constant ^{H). Similar computations can 
be attempted for the mixed derivatives, which are all equal to some common 
value rj{H) at because of G's symmetry, and we will see that the sign of 
r]{H) is irrelevant. We can now write, using Taylor's formula, 

G{x, y, z, w) = -f{H) {x^ +y'^ + z^ + 'u?) 

+ ri[H) {xy + xz + xw + yz + yw + zw) 

I i' 2 I 2 I 2 I 2\ 

+ o[x +y + z +w). 

By taking + + + vJ^ sufficiently small [this corresponds to restricting 
\i — j\ and other differences to being larger than some value m = m{H), 
whose corresponding "diagonal" terms not satisfying this restriction are 
dealt with as usual], we get, for some constant 0{H) > 0, 

G{x,y,z,w) > 9{H) (x^ + y'^ + z'^ + w'^) + rj{H) {xy + xz + xw + yz + yw + zw) . 

Let us first look at the terms in (A. 2) corresponding to x'^ + y'^ + z'^ + w'^ . 
These are collectively bounded below by the same sum restricted to i = 
j + m, which equals 

The fact that the final factor contains (i — instead of [i — which 
we had, for instance, in (A. 3) in the proof of the lemma, does not help us. 
In particular, calculations identical to those following (A. 3) show that the 
above is larger than 

2N'^"-^g{m/N), 

which does not go to if > 2/3 since g{0) calculated from (A. 4) is positive. 

For the terms in (A. 2) corresponding to xy + xz + xw + yz + yw + zw, 
considering, for instance, the term xy, similar computations to those above 
lead to the corresponding term in S being equal to 

i'=m. k'=m 



E 

{i'+0,i,k'+l,l)&D 
l<j<N-i';l<l<N-k' 



+ -l\-\j 



H-l 



" jv— i /-i — 2'/iV /-i — K /iV 

~2iV2^-2 V- V- (i/fc/)//-2 / / ^y^^ 

■I 1.1 "'0 Jo 

i—m k—m 

( i' k' \ 

X {z-w){ — + y- — -w\{y-w] 



H-l 
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which evidently tends to provided H <1. 

We conclude that if H > 2/3, then \\gN 'S'l S'Af lli2([o i]2) does not tend to 
and, by the Nualart-Ortiz-Latorre criterion [Theorem 2.1 part (iii)], U2, as 
defined in (3.32), does not converge in distribution to a normal. Hence, we 
can guarantee that, provided H > 2/3, the adjusted variation in Theorem 
3.3 does not converge to a normal. Thus, the normality of our adjusted 
estimator in Theorem 4.1 holds if and only if H £ (1/2,2/3). 
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